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Increased aridity and human population have reduced tree cover in parts of the African Sahel and 
degraded resources for local people. Yet, tree cover trends and the relative importance of climate and 
population remain unresolved. From field measurements, aerial photos, and Ikonos satellite images, we 
detected significant 1954-2002 tree density declines in the western Sahel of 18 ± 14% (P = 0.014, 
n = 204) and 17 ± 13% (P = 0.0009, n = 187). From field observations, we detected a significant 1960- 
2000 species richness decline of 21 ± 11% (P = 0.0028, n = 14) across the Sahel and a southward shift of 
the Sahel, Sudan, and Guinea zones. Multivariate analyses of climate, soil, and population showed that 
temperature most significantly (P < 0.001) explained tree cover changes. Multivariate and bivariate tests 
and field observations indicated the dominance of temperature and precipitation, supporting attribution 
of tree cover changes to climate variability. Climate change forcing of Sahel climate variability, partic- 
ularly the significant (P < 0.05) 1901-2002 temperature increases and precipitation decreases in the 
research areas, connects Sahel tree cover changes to global climate change. This suggests roles for global 
action and local adaptation to address ecological change in the Sahel. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Introduction 

In the African Sahel, people depend on trees for wood, shade, 
protection of soil fertility, and other natural resources and services. 
In the 20th century, rainfall in the Sahel declined 20—30% (Hulme 
et al„ 2001 ; Dai et al„ 2004), the most negative rainfall trend in 
the world in the instrumental record (Trenberth et al., 2007), and 
mean temperature increased up to 1.3 °C (Hulme et al„ 2001; 
Trenberth et al., 2007). As a result, aridity in the Sahel increased, 
degrading tree cover in many areas (Poupon, 1980; Gonzalez, 2001 ; 
Wezel and Lykke, 2006; Maranz, 2009) and leading to the death of 
a quarter million people in a 1968—1972 drought (U.S. Centers for 
Disease Control, 1973). This prompted countries to adopt, in 1994, 
the United Nations Convention to Combat Desertification. 

The rainfall decline and temperature increase comprise the two 
principal forms of 20th century Sahel climate variability. Increasing 
sea surface temperature governs the decline in rainfall (Giannini 
et al., 2003; Held et al., 2005; Hegerl et al., 2007; Shanahan et al., 
2009), while reduction of vegetation cover amplifies the decline 
through positive feedbacks between precipitation and vegetation 
via reduced evapotranspiration (Zeng et al„ 1999) and increased 
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albedo (Charney, 1975). Although lake sediment data indicate that 
the recent Sahel drought falls within the magnitude of West African 
droughts during the past 2500 years (Shanahan et al., 2009), 
increased atmospheric concentrations of greenhouse gases during 
the 20th century have increased sea surface temperatures, leading 
to the recent reduction in Sahel rainfall (Held et al., 2005; Zhang 
et al., 2007). Increased atmospheric concentrations of greenhouse 
gases have increased land temperatures globally, including in the 
Sahel (Hegerl et al., 2007). Thus, instrumental measurements and 
computer modeling attribute the two principal forms of 20th 
century Sahel climate variability to anthropogenic climate change 
(Wang and Eltahir, 2002; Held et al„ 2005; Biasutti and Giannini, 
2006; Hegerl et al., 2007; Zhang et al., 2007). 

In West and Central Africa, increasing rainfall toward the 
Equator (Fig. la) and decreasing temperature toward the coasts 
(Fig. lb) increase vegetation productivity (Fig. lc) and differentiate 
vegetation into three latitudinal ecological zones of increasingly 
mesic tree species and increasing biodiversity: the Sahel, the Sudan, 
and Guinea (Aubreville, 1950) (Fig. Id). The Sahel, the Sudan, and 
Guinea are biome-level zones of, respectively, savanna with shrubs 
and thorny tree species with bipinnately compound leaves, tropical 
woodland of trees with pinnately compound leaves and dry fruits, 
and tropical deciduous broadleaf closed-canopy forest of trees with 
simple broad leaves and moist fruit. In the Sahel and Sudan, indi- 
vidual woody plants of the same species grow in a continuum of 
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Fig. 1 . Environmental conditions and research areas, from spatial analyses by the authors of data from the listed sources, (a) Precipitation (1901-2002 mean) (Mitchell and Jones, 
2005) increases toward the Equator and (b) temperature (1901-2002 mean) (Mitchell and Jones, 2005) decreases toward the Atlantic Ocean, creating patterns of (c) vegetative 
productivity (1982-2008 mean NDVI) (Tucker et al., 2005) and differentiating vegetation into (d) latitudinal ecological zones (Aubreville, 1950; FAO, 2001). In the (e) aerial photo 
and Ikonos research areas (Njoobeen Mbataar [NJ], Fete Ole [FE], Wolum [WO]), (f) precipitation (Mitchell and Jones, 2005) is significantly correlated to latitude (north-south 
transect from 19°N, 15°W to 13°N, 15°W) and (g) temperature (Mitchell and Jones, 2005) is significantly correlated to longitude (east-west transect from 16 N, 12°W to 16°N, 
17° W). 


forms ranging from short multi-stemmed plants with a low crown 
(shrubs) to plants with a single trunk and elevated crown (trees). 
For brevity, trees and shrubs are hereafter referred to as trees. 

Field research across the Sahel has documented decreases in 
tree density and changes in tree species richness in the last half of 


the 20th century that were not directly caused by local people. Tree 
density declined at sites in Senegal (Poupon, 1980; Gonzalez, 2001 ; 
Vincke et al., 2010) while arid species expanded and mesic species 
declined in Burkina Faso (Rasmussen et al., 2001 ; Wezel and Lykke, 
2006; Maranz, 2009), Niger (Wezel and Lykke, 2006), and Senegal 
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(Gonzalez, 2001; Wezel and Lykke, 2006; Vincke et al., 2010). 
On the other hand, research showed no decrease in tree 
density at sites in Burkina Faso (Rasmussen et al., 2001), Mali 
(Hiernaux et al., 2009), and the Republic of the Sudan (Schlesinger 
and Gramenopoulos, 1996). Agricultural expansion and livestock 
grazing, not climate, reduced tree density at sites in Mauritania 
(Niang et al., 2008) and Senegal (Lericollais, 1989). Recent natural 
regeneration by farmers increased tree density at a site in Burkina 
Faso (Reij et al„ 2005). 

Some studies have attempted to assess Sahel vegetation cover 
trends using the normalized difference vegetation index (NDVI; 
Tucker, 1979), derived from the Advanced Very High Resolution 
Radiometer sensor of the National Oceanic and Atmospheric 
Administration polar-orbiting meteorological satellite series (e.g. 
Anyamba and Tucker, 2005; Seaquist et al., 2006). Others have used 
rain-use efficiency (RUE; Le Houerou, 1984), calculated using NDVI 
(e.g. Nicholson et al., 1998; Prince et al., 1998; Hein et al., 2011). 

Two characteristics of the NDVI data series constrain the inter- 
pretation of results of those studies. First, NDVI in the Sahel mainly 
tracks interannual variability of herbaceous vegetation. NDVI is an 
index of the fraction of photosynthetically active radiation absor- 
bed by plants, which is proportional to green leaf area. In the Sahel, 
green leaf area of the herbaceous layer, composed mainly of annual 
grasses, substantially exceeds the green leaf area of trees and 
shrubs, which allocate a higher fraction of biomass to wood, 
undetected by NDVI. Second, annual NDVI coverage began in 1982. 
A drought in the Sahel from 1982 to 1984 and a slight increase of 
rainfall in the 1980s and 1990s led to short-term increases in NDVI 
(Anyamba and Tucker, 2005; Seaquist et al., 2006) and RUE 
(Nicholson et al., 1998; Prince et al„ 1998) in some areas, showing 
possible increases in vegetation cover (Olsson et al., 2005). Because 
of the particular characteristics of the NDVI data series, however, 
NDVI can only show that herbaceous vegetation increased and only 
after 1982. 

In contrast, tree density and biodiversity indicate trends of 
vegetation cover over longer time periods than NDVI because trees 
often withstand years of stress before dying (Allen et al., 2010). 
Furthermore, changes in tree species composition can reveal 
changes in environmental factors even if total tree density does not 
change. Aerial photographs and high-resolution satellite images 
can provide data on tree densities. Structured surveys of traditional 
knowledge and field surveys can provide data on past and current 
tree species distributions. 

The decline in rainfall and rise in temperature in the Sahel 
occurred at approximately the same time as a substantial increase 
in population. Population in the Sahel grew during the 20th 
century, and, from 1960 to 1995, increased from 21 million people 
to 46 million in the nine Sahel countries of Central and West Africa 
(OECD, 1998). The increases in aridity and population have exerted 
pressure on tree cover in the Sahel, although the relative impor- 
tance of these two factors remains unresolved. Improved knowl- 
edge of vegetation change in the Sahel could increase the 
effectiveness of programs to address the vulnerability of people and 
management of natural resources in the region. We investigated 
these issues with two research objectives: (1) to detect long-term 
changes in the density and species richness of trees and shrubs in 
the Sahel and (2) to attribute causes of detected changes. 

2. Materials and methods 

2.3. Research areas for tree density 

We quantified tree density in three research areas in the 
western Sahel (Fig. le), where precipitation shows a significant 
latitudinal gradient (Fig. If) and temperature shows a significant 


longitudinal gradient (Fig. lg). Each area is a 10 km (east— west) by 
20 km (north— south) rectangle. The research areas (and center 
latitude and longitude) are: Njoobeen Mbataar, Senegal (15.80°N, 
16.25°W), Fete Ole, Senegal (16.23°N, 15.10°W), and Woium, 
Mauritania (16.25°N, 13.67°W). They are spaced at ~ 150 km along 
approximately the same latitude. Past research examined trees and 
shrubs at Njoobeen Mbataar (Gonzalez, 2001) and Fete Ole 
(Poupon, 1980). The Njoobeen Mbataar area consists of agricultural 
fields with trees, interspersed with shrubland and woodland, 
inhabited by the Wolof ethnic group. The Fete Ole area consists of 
grassland used by the Fulbe ethnic group for grazing livestock, 
interspersed with shrubland and woodland. The Woium area 
consists of grassland with single shrubs and trees, used by the Hal 
Pulaar-en ethnic group for grazing livestock. Although some 
farmers use low-intensity fires to clear fields of weeds before 
planting, wildfires are not as common as in areas to the south 
and result mainly from anthropogenic, not natural, ignition 
(Poupon, 1980). 

2.2. Research village areas for tree species richness 

We quantified tree species richness in 14 village areas in five 
countries, sampling much of the longitudinal extent of the Sahel 
(Fig. la— d). During field trips to rural areas, local staff of the Reseau 
de Systemes d’Alerte Precoce contre la Famine (Famine Early Warning 
Systems Network) selected villages not previously known to them 
or the authors. We sought rural villages that met the following 
criteria: medium population (~ 150— 300 people), agriculture or 
livestock herding as the predominant occupation, minimal urban 
infrastructure, absence of deep-bore wells that would attract large 
herds of livestock that could alter vegetation (Vincke et al., 2010), 
and presence of elder residents with detailed knowledge of natural 
history. 

2.3. Field measurements of trees 

To characterize the trees and shrubs of the three research areas 
and calibrate remote sensing data to field observations, we delin- 
eated 12 1-ha square quadrats located at the centers of the four 
quarters (NW, NE, SE, SW) of each research area. In May 2002, we 
identified the species (Table Al ) and counted all live and dead trees 
and shrubs in the quadrats. We measured all tree and shrub heights 
(Ai), diameter at breast height ( h = 1.3 m), diameter at h = 40 cm, 
diameter at ground level (shrubs and trees of h < 40 cm), and 
number of branches. In November 2005, we re-inventoried the 
Njoobeen Mbataar NW quadrat. We created geographic informa- 
tion system (CIS) point files of field tree locations. 

2.4. Tree density 

We determined current and past tree densities from Ikonos 
satellite images and historical aerial photographs, respectively, 
calibrated and validated by the field measurements of trees. The 
Ikonos satellite acquired cloud-free scenes of the Njoobeen Mba- 
taar (11:56 AM UTC, February 17, 2002), Fete Ole (11:38 AM UTC, 
February 22, 2002), and Woium (11 :37 AM UTC, February 22, 2002) 
research areas from a sun-synchronous orbital altitude of 680 km. 
We timed acquisition to match the months of historical aerial photo 
missions. February is the middle of the dry season in the Sahel, 
providing high contrast of trees over bare ground. Ikonos captured 
one band of panchromatic data at 1 m spatial resolution and four 
bands of multi-spectral data at 4 m spatial resolution, with each 
band as a GeoTiff file at 11 bit data depth. Space Imaging, Inc. 
provided geometrically corrected Ikonos data with a horizontal 
accuracy (root mean square error) of 8 m. We used the 
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panchromatic band, setting display histograms to view the images 
with approximately the same brightness and contrast observed in 
the field. 

We acquired 31 aerial photos on paper: Njoobeen Mbataar, 
February 1954, Mission AOF 083 ND-28-XX 1954 (9 photos): Fete 
Ole, ~ February— April 1954, Mission AOF 088 NE-28-1II 1954 
(8 photos); Wolum, ~ February— April 1954, Mission AOF 090 NE- 
28-V 1954 (8 photos): from the Institut Geographique National 
(IGN), France at 1: 50,000 scale; and Njoobeen Mbataar, March 
1989, Japan International Cooperation Agency (JICA) Mission 1989 
(6 photos), from the Service Geographique National, Senegal at 1: 
60,000 scale. We scanned the IGN photos at 500 dots otT 1 and the 
JICA photos at 600 dots cm -1 to produce 8 bit tiff images of 1 m 
spatial resolution. We adjusted display histograms to the same 
brightness and contrast as the Ikonos images. We geographically 
registered the aerial photographs to the Ikonos images using 5937 
control points and formed mosaics that matched the extent and 
1 m spatial resolution of each Ikonos image. 

Using ENV1 4.6.1 CIS software, we viewed all images at the same 
brightness on a MacBook Pro color liquid crystal display with 
dimensions 33 cm by 20.7 cm ( 1440 pixels by 900 pixels) and 32 bit 
color depth. Using the Ikonos image centers as base points, we 
divided each research area into 1 km 2 cells and delimited a 1 ha 
square sample area at the center of each cell. The Njoobeen Mbataar 
area had 200 sample areas, the Fete Ole area had 182 sample areas 
(because of incomplete aerial photo coverage on the east edge), and 
the Wolum area had 198 sample areas (because of incomplete 
aerial photo coverage in the southwest corner and because one 
sample fell inside a village). Each area also had four field quadrats. 
To confirm the independence of the samples, we calculated global 
Moran’s / (Anselin, 1995), an index of spatial autocorrelation, on 
tree density for each research area and year. 

To calibrate aerial photo and Ikonos image data to the field 
measurements of trees, we established the image properties of field 
trees of It > 3 m. For each aerial photo mosaic and Ikonos image, we 
determined thresholds of minimum and mean tree center dark- 
ness, crown area, and contrast between the crown and the 
surrounding area of trees of h > 3 m in the 12 field quadrats and, to 
generate an adequate sample size for the Wolum area, six adjacent 
1-ha sample areas in the Wolum 2002 Ikonos scene. 

To complete the calibration of aerial photo and Ikonos image 
data to the field measurements of trees, we compared the number 
of trees appearing in Ikonos images viewed at 1:7500 scale to the 
trees of /i > 3 m measured in the 12 field quadrats. Previous vali- 
dation of tree counts derived from 1989 aerial photos against field 
inventories showed no statistical difference between dots visible 
at 1:7500 scale and trees of h > 3 m measured in the field 
(Gonzalez, 2001). 

At 1:7500 scale, we viewed each sample area and quadrat in 
each aerial photo mosaic and Ikonos image. Using the aerial photo 
and Ikonos thresholds (Table A2) to determine objects to include as 
trees, we counted trees of h > 3 m and created CIS point files of 
their locations. To reduce viewing errors, we repeated the tree 
counts three times for each research area and year. 

We calculated 2002 tree densities by dividing raw tree counts 
from Ikonos by the fraction of field trees visible in Ikonos images 
(Table A3). To calculate tree mortality, we tracked individual trees 
through all time periods. To calculate significance of net differences 
in tree density, we used heteroscedastic t-tests on paired data. 

2.5. Tree species richness 

From January 2000 to May 2002, we visited the 14 research 
villages. Each village identified male and female elders ~60 years 
old whom we interviewed as a group. We proceeded through a list 


of 161 Sahel, Sudan, and Guinea tree and shrub species present in 
the Sahel (Table Al ), identified by local name, and asked them to 
identify species as present or absent in the village lands in 1960, 
a year memorable as the year of independence of each country from 
France, and in 2000. We used botanical characteristics to verify the 
identity of species. We also asked them to identify any tree or shrub 
species that we had not already listed. Before and after speaking 
with the villagers, we independently surveyed village lands to 
validate their information on species distributions in 2000. 

We calculated species richness (number of species) and tree 
species composition (fraction of all species that were characteristic 
of each ecological zone). To normalize results among village lands 
of a range of surface areas, we calculated 1960—2000 changes as 
a fraction of 1960 values. To calculate significance of differences, we 
used heteroscedastic t-tests on paired data. 

2.6. Attribution of causes of change 

To attribute causes of changes in tree density and tree species 
richness, we conducted four canonical correlations analyses 
(CCA; Hotelling, 1936), including one for each research area and one 
for the group of research village areas. We used CCA to examine the 
relative influence of four fundamental factors — temperature, 
precipitation, soil, and population— on tree density and tree species 
richness. Previous research has applied CCA to the attribution of 
causes of changes in Sahel tree density and species composition 
(Gonzalez, 2001) and to other direct comparisons of factors influ- 
encing the distribution of tree species (e.g. Carleton, 1984; Pelissier 
et al„ 2002). To assess the robustness of CCA results and further 
examine relationships between dependent and independent vari- 
ables, we also calculated bivariate correlations for each pair of 
variables and the significance of the correlations. 

For Njoobeen Mbataar, we conducted CCA on nine variables: 
tree density in 1954, 1989, and 2002, tree density change as 
a fraction of original tree density for the periods 1954—1989, 
1989—2002, and 1954—2002 (dependent variables); 1950—2000 
mean temperature and precipitation (Hijmans et al., 2005), pop- 
ulation density in 2008 (Dobson et al., 2000) (independent vari- 
ables). Only one soil type occurs in the Njoobeen Mbataar area (FAO 
et al., 2009), reflecting low variation in the area’s soil characteris- 
tics, so soil was not included in the Njoobeen Mbataar CCA. For Fete 
Ole and Wolum, we conducted CCA on seven variables: tree density 
in 1954 and 2002, tree density change from 1954 to 2002 as 
a fraction of 1954 tree density (dependent variables); 1950—2000 
mean temperature and precipitation (Hijmans et al„ 2005), topsoil 
(0—30 cm depth) organic carbon fraction (FAO et al., 2009), pop- 
ulation density in 2008 (Dobson et al., 2000) (independent vari- 
ables). We projected all data from the geographic reference system 
at a spatial resolution of 30 s to the equal-area Universal Transverse 
Mercator (UTM) projection at a spatial resolution of 1 km. 

For the set of Sahel village areas, we conducted CCA on seven 
variables: tree species richness in 1960 and 2000, change in tree 
species richness as a fraction of 1960 richness (dependent vari- 
ables); 1950—2000 mean temperature and precipitation (Hijmans 
et al., 2005), topsoil (0—30 cm depth) organic carbon fraction 
(FAO et al., 2009), population density in 2008 (Dobson et al., 2000) 
(independent variables). We used the UTM data at 1 km spatial 
resolution. To estimate parameter values for each terroir villageois 
(each set of village lands), we calculated average values for the 
4 km 2 (four pixels) around each village, the approximate extent 
estimated by village residents and an area consistent with previous 
analyses (Gonzalez, 2001 ). 

To examine if tree density and species richness change were 
attributable to climate change, we determined if each research area 
and research village was located in an area of 20th century climate 
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change. This contributed to a method of joint attribution previously 
used to connect historical ecological changes to anthropogenic 
climate change (Rosenzweig et al., 2008). We used the University of 
East Anglia Climate Research Unit TS 2.1 data set of 1901—2002 
climate (Mitchell and Jones, 2005), from which the Intergovern- 
mental Panel on Climate Change has detected significant changes in 
temperature and precipitation (Trenberth et al„ 2007) and attrib- 
uted the detected changes to anthropogenic climate forcing (Hegerl 
et al., 2007). We projected the data from the geographic reference 
system at a spatial resolution of 0.5° to Lambert Azimuthal Equal- 
Area projection at a spatial resolution of 50 km and derived 
climate trends by linear least squares regression of observed 
mean annual temperature and annual precipitation versus time 
(Gonzalez et al., 2010). 

2.7. Limitations of the methods 

Aerial photos and Ikonos satellite images differed in quality. 
Ikonos images showed finer detail, demonstrated by smaller 
thresholds for detecting trees (Table A2). This difference would 
tend to cause undercounting of trees in the aerial photos, implying 
that tree density decreased even more than estimated. Calibration 
of Ikonos data to field observations showed that tree counts derived 
from Ikonos underestimated field densities (Table A3). Conse- 
quently, we calculated 2002 tree densities by dividing raw tree 
counts by the fraction visible in each area, further narrowing the 
difference between tree densities derived from aerial photos and 
Ikonos images. We validated our tree counts at Fete Ole against 
field inventories of the same area by Poupon (1980). 

For the species surveys, recollections of the presence or absence 
of species served as a proxy for non-existent data on historical 
distributions. Although recollections are inexact, we systematically 
recorded just one piece of binary information (presence = 1, 
absence = 0) for each species for each year and used a past year 
marked by a memorable event. We verified the accuracy of 
villagers’ responses on current species through our own field 
observations. Previous validation of local knowledge of tree species 
gave an error rate of 1% of data points (Gonzalez, 2001), equivalent 
to approximately one species for each of the 14 Sahel village areas. 

For the attribution of causes, we recognize that our choice of 
four independent variables is a simplification. We sought to 
account for the myriad environmental and socio-economic vari- 
ables involved by selecting the fundamental factors most influ- 
encing tree cover in the Sahel. Availability of spatial data also 
constrained our choice of variables. Furthermore, we have sought 
to account for complex interactions among variables by using CCA, 
a multivariate statistical method that systematically examines the 
relationships of every pair of variables and calculates the influence 
of each independent variable on the observed variability of all the 
dependent variables. Although CCA provides a comprehensive 
comparison of the relationships of sets of dependent and inde- 
pendent variables, canonical loadings prove correlation not 
causation. 

Although data on mean climate are available at 1 km spatial 
resolution (Hijmans et al., 2005), data on climate trends are only 


available at 50 km spatial resolution (Mitchell and Jones, 2005; 
Gonzalez et al., 2010). Although data on population density are 
available (Dobson et al., 2000), consistent data on population 
trends are not available. In addition, the LandScan population 
density data (Dobson et al., 2000) were only available for 2008 
because the developers do not distribute data from previous years. 

3. Results 

3.3. Tree density 

Field counting and measurement of 1503 live trees in the 12 
quadrats (Table 1) showed that Fete Ole had the highest tree 
density while Njoobeen Mbataar had the highest species richness 
and highest fractions of trees of Sudan and Guinean species. 
Because we only found four dead trees and because all of them 
were small (h < 40 cm), we analyzed data from only the live trees. 
Validation of our tree counts against a field inventory of the 25 ha 
site at Fete Ole (Poupon, 1980) shows general agreement, with 
a 29% decrease in tree density (1972—1976) derived from field data 
and a 20% decrease (1954—2002) derived from aerial photos and 
Ikonos. Spatial autocorrelation was low in all research areas, with 
the following values of global Moran’s 3: 0.163, 0.067, 0.016 (Njoo- 
been Mbataar 1954, 1989, 2002); 0.145, 0.052 (Fete Ole 1954, 2002); 
0.065, 0.035 (Wolum 1954, 2002). 

We found significant decreases in tree density at Njoobeen 
Mbataar from 1954 to 1989 and from 1954 to 2002 and at Fete Ole 
from 1954 to 2002 (Fig. 2). At both Njoobeen Mbataar and Fete Ole, 
tree density fell one-fifth from 1954 to 2002 and fell in over half the 
samples (Table A4). At Wolum, tree density fell in one-third of the 
samples, but net change in average density was not significant. By 
tracking individual trees (Fig. 3), we found that tree mortality rates 
ranged from 0.11 ± 0.01 decade -1 at Wolum to 0.33 ± 0.04 
decade -1 at Njoobeen Mbataar (Table A4). Tree density and tree 
density change displayed substantial spatial variation across each 
research area (Figs. Al— A3). 

3.2. Tree species richness 

Field surveys recorded 116 tree and shrub species in the 14 
village areas (Table Al), with the following species composition: 
Sahel species 29%, Sudan species 39%, Guinean species 32%. The 
family Mimosaceae had the highest number of species of each 
ecological zone, with Capparaceae (Sahel), Caesalpiniaceae (Sudan), 
and Moraceae (Guinea) representing the second highest number of 
species of each ecological zone. Tree species richness declined in 12 
village areas and remained unchanged in two village areas 
(Table A5). 

We found a significant decrease in tree species richness in the 
group of village areas from 1960 to 2000 (Fig. 2). Tree species 
richness of each of the three ecological zones fell significantly, with 
the rate of decrease in Guinean species double the rate of decrease 
in Sahel species (Table 2). As a fraction of all tree species, Sahel 
species increased significantly while Guinean species decreased 
significantly (Table 2). 


Table 1 

Tree densities and tree species (mean ± standard deviation [SD]) by research area, listed from west to east. 


Research area 

Country 

Tree density (trees ha 1 ) 



Species richness 
(species) 

Species composition (fraction of all trees) 

Trees 

counted 

n 

h < 40 cm 

40 cm < h < 3 m 

h > 3 m 

All trees 

Mean 

Total 

Sahel 

Sudan 

Guinea 

Njoobeen Mbataar 

Senegal 

33 ± 18 

87 ± 24 

9.5 ± 8.8 

130 ± 33 

10 ±2 

18 

0.55 ± 0.20 

0.40 ±0.18 

0.05 ± 0.02 

518 

4 

Fete Ole 

Senegal 

12 ± 7.2 

134 ± 77 

21 ± 11 

167 ± 89 

6±3 

13 

0.99 ± 0.30 

0.01 ± 0.03 

0 

669 

4 

Wolum 

Mauritania 

12 ± 9.0 

67 ± 69 

0.2 ± 0.5 

79 ± 79 

4±3 

9 

0.83 ± 0.34 

0.17 ± 0.34 

0 

316 

4 
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Fig. 2. Changes in tree cover and climate, (a) Decrease in tree species richness of 
village areas, (b) Decrease in tree density at Fete Ole. (c) Changes in tree density at 
Njoobeen Mbataar. (d) Changes in tree density at Wolum. (e) Mean temperature 
(calculated from data from Mitchell and Jones (2005); thin line) from 1901 to 2002 at 
St. Louis, Senegal (16.1°N, 16.5°W; location shown in Fig. la), running five-year mean 
(thick line), and significant trend (dotted line) from least squares linear regression, 
(f) Annual precipitation (thin line; calculated from data from Aubreville (1938) and 
from National Oceanic and Atmospheric Administration, Global Historical Climatology 
Network, http://www.ncdc.noaa.gov/oa/climate/ghcn-monthly/index.php) from 1854 
to 2008 at St. Louis, Senegal, running five-year mean (thick line), and significant trend 
(dotted line) from least squares linear regression. 

3.3. Attribution of causes of change 

Temperature had the highest rotated loading for canonical 
function (CF) 1 for all four CCAs and precipitation had the second 
highest rotated loading for CF 1 for three CCAs (Table A6). All cor- 
responding loadings were significant (P < 0.0001 ). For CF 2, topsoil 
organic carbon had the highest rotated loading for Wolum and the 
Sahel villages, population had the highest rotated loading for 
Njoobeen Mbataar, and precipitation had the highest rotated 
loading for Fete Ole (Table A6). 


In bivariate correlation tests, tree density change (1954—2002) 
at Njoobeen Mbataar showed significant correlation with temper- 
ature (r 2 = 0.21, P = 0.003, n = 200) and precipitation (r 2 = -0.21, 
P = 0.002, n = 200). Species richness change (1960—2000) in the 
Sahel villages showed significant correlation with temperature 
(r 2 = 0.67, P = 0.009, n = 14). Tree density and tree species richness 
changes did not show significant bivariate correlations with topsoil 
organic carbon or population density. 

According to government censuses, population grew in the 
administrative regions of all three research areas: Njoobeen Mba- 
taar, 0.02 year -1 (Region de Louga, 1976—2002, http://www.ansd. 
sn); Fete Ole, 0.01 year -1 (Region de St. Louis, 1976—2002, http:// 
www.ansd.sn); Wolum, 0.03 year -1 (Moughataa de M’Bagne, 
1977—2000, http://www.ons.mr). Population growth was highest 
in Wolum, the only area where tree density did not significantly 
change. 

The detected changes in tree density and species richness 
occurred in areas of climate change (Table 3, Fig. 4). At all three 
research areas, temperature increased and rainfall decreased 
significantly in the 20th century. The meteorological station at St. 
Louis, Senegal, 30—300 km west of the research areas (location 
shown in Fig. la), has recorded the longest rainfall data series in 
Africa and documented the rainfall decline, showing a significant 
linear trend of —1.4 mm year -2 from 1854 to 2008 (Fig. 2). 
Temperature in the St. Louis region showed a significant linear 
trend of 0.7 °C century -1 from 1901 to 2002 (Fig. 2). Twentieth 
century temperature increased in 13 village areas, with significant 
increases in 10 village areas (Table 3, Fig. 4). Rainfall decreased in 12 
village areas, with significant decreases in three village areas. 
Whereas temperature increased and precipitation decreased for 
most of the areas, covariance of temperature change with precipi- 
tation change was low across the Sahel (covariance = -0.02 for the 
land area of Fig. 4). 

In summary, during the 20th century, average tree species 
richness decreased significantly in the research village areas 
(Fig. 2a), tree density decreased significantly in two research areas 
(Fig. 2b and c), but remained stable in the most arid area (Fig. 2d), 
temperature increased significantly (Fig. 2e), and precipitation 
decreased significantly (Fig. 2f). 

4. Discussion 

4.1. Significant declines in tree density and species richness 

The significant changes in tree density at Njoobeen Mbataar and 
Fete Ole provide evidence of tree density declines in the western 
Sahel. The decreases in tree density occurred at rates close to the 
rates previously found in parts of the two areas (Gonzalez, 2001; 
Vincke et al., 2010). Our sampling was more intensive in the 
Njoobeen Mbataar area and more extensive in the Fete Ole area 
than previous research (Poupon, 1980; Gonzalez, 2001; Vincke 
et al., 2010). Although the change at Wolum was not significant, 
tree density at Wolum was one-fifth and one-tenth the densities at 
Njoobeen Mbataar and Fete Ole, respectively, and 1950—2000 
mean temperature was 2 °C and 4 °C warmer, respectively. It is 
possible that arid conditions had already thinned trees to low 
densities before 1954. 

A significant decline in tree density at Njoobeen Mbataar from 
1954 to 1989, followed by no significant change from 1989 to 2002 
(Fig. 2c), is consistent with observations by local residents and data 
from Fete Ole (Poupon, 1980) of substantial tree mortality in the 
1970s and 1980s. The 154-year record of precipitation at St. Louis, 
Senegal (Fig. 2f) shows a significant rainfall decline that hit its nadir 
in the period 1968—1973, when rainfall fell to approximately one 
standard deviation below the long-term average five out of six 
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Fig. 3. Aerial photos and Ikonos panchromatic images of a 1 km 2 section of the Njoobeen Mbataar area and trees identified in the 1 ha sample area around 15.75°N, 16.25°W. 


years. This is a possible threshold of resilience to drought stress for 
Sudan and Guinean tree species. 

The significant declines in tree density in Senegal (Poupon, 
1980; Gonzalez, 2001; Vincke et al., 2010; this research) contrast 
with cases of no significant change or increases in tree or shrub 
density at sites in Burkina Faso, Mali, and the Republic of the Sudan 
(Schlesinger and Gramenopoulos, 1996; Rasmussen et al., 2001; 
Hiernaux et al., 2009). This suggests that the heterogeneity of 
environmental conditions in the Sahel has led to uneven changes in 
tree cover. In addition, the significant decline of tree density at 
Njoobeen Mbataar and Fete Ole from 1954 to 2002 contrasts with 
the lack of significant change in NDV1 in the two areas from 1982 to 
1999 (Seaquist et al., 2006). This demonstrates that tree density can 
reveal long-term trends undetectable by NDV1 because of the 
starting year of the NDVI time series and its closer relationship, in 
the Sahel, to herbaceous cover rather than tree cover. 

The significant decrease in tree species richness in the group of 
14 village areas and the replacement of mesic species by xeric 
species is consistent with changes in tree species at numerous sites 
in Burkina Faso, Niger, and Senegal (Gonzalez, 2001; Rasmussen 
et al., 2001; Wezel and Lykke, 2006; Maranz, 2009). The two 
village areas where we found no change in tree species richness, 
Fabugu and Kaylaroom, are located in the floodplain of the Niger 
River and the basin of Lake Chad, respectively, suggesting that 


a high water table maintained favorable soil moisture, conserving 
mesic species. 

In the remaining research villages, the decrease in Guinean 
species as a fraction of all species and the increase in Sahel species 
as a fraction of all species indicate a shift of the Sahel zone south- 
ward into areas of increased aridity and a concomitant shift of the 
Guinean zone southward toward more humid areas, a vegetation 
shift previously documented in Senegal (Gonzalez, 2001). The 
results indicate a retraction of Guinean species and southward shift 
of the Sahel-Sudan and Sudan-Guinea ecotones. 

In summary, the decreases in tree density in the western Sahel 
and tree species richness across the Sahel together provide 
evidence of a decline of tree cover in the Sahel in the last half of the 
20th century. Two different types of tree data — density and species 
richness— drawn from different samples both showed significant 
declines. Tree density at Njoobeen Mbataar and Fete Ole and tree 
species richness in the research villages declined at similar rates. 

4.2. Attribution to climate variability and joint attribution 
to climate change 

Canonical correlations analyses indicate that temperature vari- 
ability dominated all factors in explaining observed patterns of tree 
density and species richness. Precipitation variability was almost as 


Table 2 

Tree species richness, tree species composition, and changes (mean ± SD). 


Ecological 

zone 

Species richness 

Species richness change 



Species composition 

Species composition change 

n 

(species) 

1960 

2000 

(species) 

(fraction) 

(decade ’) 

P 

(fraction of all species) 
1960 2000 

(fraction of 
all species) 

P 

Sahel 

17 ± 2 

15 ±2 

-2 ± 2 

-0.13 ±0.10 

-0.03 ± 0.02 

0.0140 

0.38 ± 0.07 

0.43 ± 0.07 

0.04 ± 0.03 

0.0047 

14 

Sudan 

18 ±3 

15 ±3 

-3 ± 3 

-0.18 ±0.12 

-0.04 ± 0.03 

0.0188 

0.39 ± 0.03 

0.40 ± 0.03 

0.01 ± 0.03 

0.3243 

14 

Guinea 

12 ±4 

7 ±3 

-4 ±2 

-0.35 ±0.17 

-0.09 ± 0.04 

0.0010 

0.23 ± 0.06 

0.17 ± 0.06 

-0.06 ± 0.03 

0.0011 

14 

All 

47 ± 7 

37 ± 7 

-10 ± 6 

-0.21 ±0.11 

-0.05 ± 0.03 

0.0028 

1.00 

1.00 

- 

- 

14 
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Table 3 

Climate, soil, and population characteristics of research locations, listed from west to east. Mean (±SD) calculated for each research area (n = 200 pixels, 1 pixel = 1 km 2 ) and 
each research village area (n = 4 pixels, 1 pixel = 1 km 2 ) from data on mean climate (Hijmans et al., 2005), topsoil (depth 0-30 cm) organic carbon (FAO et al., 2009), and 
population density (Dobson etal., 2000). Climate trends (Mitchell and Jones, 2005; Gonzalez etal., 2010) given for the 50 km x 50 km pixel surrounding each research location. 
Significance of climate trends from linear regression: *P < 0.05, **P < 0.01, ***P < 0.001. 



Country 

Mean temperature 

Temperature trend 

Mean precipitation 

Precipitation 

Soil organic 

Population density 



(”C) 

(°C century -1 ) 

(mm year -1 ) 

trend (century -1 ) 

carbon (%) 

(people km 2 ) 



1950-2000 

1901-2002 

1950-2000 

1901-2002 


2008 

Research areas 
Njoobeen Mbataar 

Senegal 

25.3 ± 0.08 

0.7 *** 

330 ± 9 

-0.48*** 

0.3 

79 ± 170 

Fete Ole 

Senegal 

27.6 ± 0.03 

0.8*" 

290 ±9 

-0.42*** 

0.2 ± 0.05 

13 ±4 

Wolum 

Mauritania 

29.5 ± 0.07 

0.7*** 

260 ± 7 

-0.31* 

0.9 ± 0.3 

39 ± 110 

Research village areas 
Juude Waalo 

Mauritania 

29.3 ± 0.02 

0.7*** 

240 ± 1 

-0.31* 

0.6 ± 0.7 

480 ± 940 

Aten 

Mauritania 

29.5 ± 0.04 

0.7*** 

240 ± 1 

-0.31* 

1.2 

270 ± 460 

Fabugu 

Mali 

27.2 ± 0.03 

0.8*** 

590 ±2 

-0.24** 

0.7 

57 ±30 

Nampabuum 

Burkina Faso 

28.3 ± 0.01 

0.6*** 

600 ± 1 

-0.12 

0.8 

410 ± 630 

Banizumbi 

Niger 

29.2 ± 0.01 

0.4* 

380 ± 1 

-0.02 

0.5 

56 ± 91 

Tamaka 

Niger 

28.7 ± 0.02 

0.3 

440 ±3 

0.02 

0.4 

86 ± 53 

Dan Tsuntsu 

Niger 

27.4 ± 0.01 

<0.1 

360 ± 1 

-0.07 

0.4 

210 ±330 

Guidimouni 

Niger 

27.6 ± 0.25 

-0.01* 

320 ± 2 

-0.02 

0.5 

54 ±44 

Kaylaroom 

Chad 

27.8 ± 0.02 

0.1 

290 ± 1 

0.09 

0.4 

1 ± 1 

Buurtey Ganuun 

Chad 

29.3 ± 0.04 

0.4* 

350 ± 1 

-0.10 

1.4 

4 

Marchuut 

Chad 

28.8 ± 0.02 

0.4* 

440 ± 2 

-0.10 

1.4 

11 ±9 

Akar 

Chad 

28.8 ± 0.05 

0.4* 

380 ± 2 

-0.10 

1.4 

39 ± 65 

Kardofal 

Chad 

27.9 ± 0.03 

0.5* 

480 ± 2 

-0.15 

0.4 

4 ± 2 

Ningelin 

Chad 

28.1 ± 0.03 

o 

In 

450 ± 2 

-0.15 

1.1 

3 ±2 


important as temperature variability, based on the loadings for CF 1. 
Whereas tree density change showed significant correlation with 
temperature and precipitation at Njoobeen Mbataar and tree 
species richness change showed significant correlation with 
temperature in the research villages, tree density change and tree 
species richness change showed no correlation with topsoil organic 
carbon or population. 

Field observations and interpretation of aerial photos and Iko- 
nos images confirm that tree density is not consistently higher or 
lower close to population centers (villages or large towns). Local 
people generally do not fell live mature trees for firewood, but 


periodically coppice (cut near the root crown) young shrubs, which 
resprout for future firewood harvesting. Tree density and tree 
species richness have declined in uninhabited areas to the west of 
Njoobeen Mbataar (Gonzalez, 2001 ), further suggesting a dominant 
cause other than local human population. In addition, even though 
the region around Wolum experienced the highest population 
growth rate of the three research areas, Wolum was the only 
research area where tree density did not significantly change. 

Precipitation dominated population in explaining tree cover 
change even though the rate of population growth of the three 
research areas exceeded the rate of precipitation decline by factors 



a Precipitation change 


500 km 


tania 


Burkina Faso 


1 0N 


» Ikonos * Village 


iKonos village 

[ -0.50 +0.50 century 




- 0.2 


+1.3° C centur 


Fig. 4. Climate changes in the period 1901-2002 (Gonzalez et al., 2010) attributable to anthropogenic climate forcing (Hegerl et al., 2007). (a) Annual precipitation linear trend 
(-0.50 century -1 = 50% decrease in a century), (b) Mean temperature linear trend. 
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of five and three at Njoobeen Mbataar and Fete Ole, respectively. In 
contrast, tree density (Table A4) and rainfall (Table 3) declined at 
similar rates at Njoobeen Mbataar and Fete Ole. 

Because the CCA results, bivariate correlations, and field obser- 
vations consistently indicate the dominant influence of tempera- 
ture and precipitation variability, we attribute the detected changes 
in tree density and tree species richness to climate variability. 
A convergence of 20th century trends of increasing temperature 
and decreasing precipitation, tree density, and species richness 
(Fig. 2) illustrates the effect of climate. While our analyses indicate 
that climate factors predominate over population in explaining 
overall patterns of tree cover change across the region, we recog- 
nize that people can still affect tree density and tree species rich- 
ness at specific locations. Still, the predominance that we find of 
climate over population and soil in explaining overall patterns of 
tree cover in the Sahel is consistent with continental-scale analyses 
of the factors governing tree cover in Africa (Sankaran et al., 2005). 

We connect Sahel tree cover changes to global climate change in 
a two-step method known as joint attribution (Rosenzweig et al., 
2008): (1) attribution of tree cover changes to Sahel climate vari- 
ability (20th century temperature and precipitation change) and (2) 
attribution of Sahel climate variability to global climate change. Our 
CCA results, bivariate correlations, and field observations provide 
evidence for step 1. Our detection of changes in tree density and 
species richness in areas of 20th century climate change (Table 3, 
Fig. 4) and previous research that attributes the two principal forms 
of 20th century Sahel climate variability to global climate change 
(Wang and Eltahir, 2002; Held et al., 2005; Biasutti and Giannini, 
2006; Hegerl et al., 2007; Zhang et al., 2007) provide evidence for 
step 2. 

The combination of warmer temperatures and significant 
decreases in precipitation in the three research areas and three 
village areas increased aridity and water stress on trees. The 
detected changes in tree cover were outside the range of spatial 
variation in each area and consistent with process-level under- 
standing of the vulnerability of mesic tree species to aridity. 
Drought can kill trees through cavitation of water columns within 
the xylem and through carbohydrate deficits that reduce resistance 
to pests (Allen et al., 2010). 

The joint attribution of Sahel tree cover change to global climate 
change is consistent with cases of drought-induced tree dieback 
around the world associated with climate change (Allen et al., 
2010). Furthermore, the detected southward shift of the Sahel, 
Sudan, and Guinean zones is similar to vegetation shifts attributed 
to climate change at numerous sites around the world (Gonzalez 
et al., 2010). These results are consistent with the joint attribu- 
tion of >90% of detected ecological changes in 29,000 data series 
around the world to anthropogenic climate change (Rosenzweig 
et al., 2008). 

5. Conclusions 

Using field measurements of trees, field surveys of tree species, 
aerial photos, Ikonos satellite images, multivariate analyses of 
temperature, precipitation, soil, and population, and spatial anal- 
yses of climate change data, we have shown that: 

1. Tree density declined significantly from 1954 to 2002 in the 
western Sahel at Njoobeen Mbataar and Fete Ole, Senegal. 

2. Tree species richness declined significantly from 1960 to 2000 
across the Sahel in a sample of 14 village areas in Mauritania, 
Mali, Burkina Faso, Niger, and Chad. 

3. The detected tree cover changes are attributable to two forms of 
climate variability — increasing temperature and decreasing 
precipitation — and jointly attributable to global climate change. 


This detection and attribution of tree cover changes can help 
guide future adaptation measures for natural resource manage- 
ment in the Sahel. One adaptation measure to maintain and 
augment tree cover is natural regeneration, a traditional practice in 
which farmers and herders select small field trees, protect them, 
prune them, and raise them to maturity. In parts of the Sahel, 
natural regeneration has doubled the density of mature trees 
(Lericollais, 1989). 

Our results suggest that natural regeneration should concen- 
trate on a mix of tree species that are more drought-tolerant than 
those that might currently occupy a field. Drought-tolerant trees 
include Sahel species like Acacia raddiana and Balanites aegyptiaca, 
which, although not as highly desired by local people for fruit or 
wood as Sudan species like Mangifera indica or Guinean species like 
Bombax costatum, are more likely to survive arid conditions. In 
addition to natural regeneration, field water harvesting and other 
practices tested across the Sahel (Larwanou and Saadou, 2011 ) offer 
alternatives for future adaptation. 

Because the reduction in tree cover in the Sahel may contribute 
to wind erosion, reduced soil organic matter, and other forms of 
land degradation, it shows how climate change may exacerbate 
desertification. Consequently, adaptation measures to increase tree 
cover could help to reduce the vulnerability of local people to both 
climate change and desertification. 

Finally, the connection of Sahel vegetation change to global 
climate change suggests that global greenhouse gas emissions 
reductions should accompany any local adaptation measures in the 
Sahel. 
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